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ABSTRACT 

We co mpute the cluster auto -correlation function S,cc{f') of an X-ray flux limited sample of Abell clusters 
(XBACs, Ebeling et al. 1996). For the total XBACs sample we find a power-law fit ^cc = (r/rQ)'^ with 
To — 21.1 Mpc h^^and 7 = —1.9 consistent with the results of i? > 1 Abell clusters. We also analyze 
iccif) for subsamples defined by different X-ray luminosity thresholds where we find a weak tendency of 
larger values of rg with increasing X-ray luminosity although with a low statistical significance. In the 
different subsamples analyzed we find 21 < tq < 35 Mpc h~^ and —1.9 < 7 < —1.6. Our analysis suggests 
that cluster X-ray luminosities may be used for a reliable confrontation of cluster spatial distribution 
properties in models and observations. 
Subject headings: cosmology-clusters-correlation function 



> 

(N 

O 

00 . 

Oh 
I 

o 



1. INTRODUCTION 

Different authors have analyzed the cluster-cluster spa- 
tial two-point correlation function finding po wer-law fits 
of the form f.^fr) = (r/mV with t - -1.8 ( |Bahcall & 
Soncira 1983iriPeacock fc West 199^ ). The value of the 



cluster-cluster correlation length rg is controversial, as also 
is its dependence on cluster mass. This has been achieved 
by studying samples selected by cluster richness and the 
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-^*r — where 



naaociatcd mean inter cluatcr acparation dc — 
Uc is the m ean number density of c lusters 
West 1992| and p3ahcall & Cen 1992 argue for a univer- 



Bahcall & 



sal scaling relation for the two-point correlation function 
of rich clusters where the cluster correlation length satis- 
fies ro = OAdc- At low v alues of dc — 30 — 50 h~^ Mpc 



the APM Cluster Survey (|Dalton ct al. 1994|), the Ed 



inburgh/Durham Cluster Catalog Lumsdcn ct al 
(hereafter EDC C), and the [AbeU 1958| and [Abell , (7orwi 
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& Olowin 1989 cluster samples give similar results consis 
tent with tq = 15 — 2Qh~^ Mpc. At larger dc, however, the 
analyses rely only on the Abell catalog and o n a cluster 



samp le selected from the APM Galaxy Survey ([Croft et al 



1997 ). The results of this high richness APM cluster sam- 
ple are not consistent with th e universal scaling rel ation 



derived from Abell clusters by Bahcall & West 1992 since 
only a weak dependence of ro on dc is found in the rich 
APM cluster sample. A partial explanation for the differ- 
ent results between rich Abell and APM clusters could rely 
on the fact that the Abell catalog is subject to visual arti- 
ficial inhomogeneities in contrast wit h the automated and 
well controlled APM cluster catalog ( [Croft et al. 1997| ). It 
should be noted, however, that given the steeper correla- 
tions for richer clusters these results are not inconsistent 
with the universal relation in terms of correlation ampli- 
tude. 

T he problems of projection effects in cluster selection 
(see van Haarlem et al. 1997| ) may be strongly overcome 
by selecting clusters in the X-ray rather than the opti- 
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cal. Moreover, given the good correlation between X-ray 
luminosity and clust er mass {Lx cx: M a) found in both 
analytical so lutions (Bertschinger 1985[ ) and in numerical 
simulations ( [Navarro et al. 1995 ), an X-ray selected sam- 
ple is suitable to study the dependence of cluster spatial 
correlations on mass. In this work we explore the values 
of ro in subsamples taken from the X-ray br ightest Abell- 



type clusters of galaxies, hereafter XBACs (Ebeling et al 
1996). This sample of clusters is complete in X-ray fiux 



and we have selected subsets with different cuts in X-ray 
luminosity L^- 

2. DATA AND ANALYSIS 

The X-ray-brightest Abell- type clusters of galaxies sur- 
vey (hereafter XBACs, Ebeling et al. 1996 ) comprises 277 
objects and is a 95 % complete flux limited sample. We 
have restricted this catalog to galactic latitudes |6| > 25° 
and X-ray flux / > fcut — 5 x 10~^^ erg cm~^ s~^ in 
the 0.1 — 2.4 keV band comprising a final sample of 248 
clusters. This sample although optically selected is con- 
firmed by the X-ray emission of the intra-cluster gas thus 
excluding spurious Abell cluste rs generated by projec tion 
effects. Also, as discussed by Ebeling et al. 1996, the 



XBACs sample is unaffected to first order by the incom- 
pleteness in volume of the Abell catalog at large distances 
since missing Abell clusters of low richness would not be 
included in XBACs due to their low X-ray luminosity. 
In figure 1 we plot the X-ray luminosity of the clusters 
as a function of redshift z taken from Table 3 of EbeL 



L 



ing et al. 1996 . The smooth curve corresponds to the 
luminosity of an object with fiux fcut at redshift 2; in a 
flat cosmology. Cluster distances d we re derived using the 
standard relation (e.g., Bandage 1961 ) 



_ c(goz + (go - 1)((1 + 2goz)V2 - 1)) 
where z is the cluster redshift, ho is the Hubble constant 



2 



in units of 100 knis^^ Mpc^^ , and c is the speed of light. 
Throughout this paper we have adopted a deceleration pa- 
rameter qo = 0.5. 




Figure 1. X-ray luminosity vs redshift z for the total 
sample of of 248 clusters analyzed. The smooth curve displays 
the luminosity corresponding to flux fcm at redshift z in a flat 
cosmology. 

We calculate the cluster-cluster two-point spatial cor- 
relation function ^cc('r) cross correlating the data with a 
random catalog constructed by randomizing the angular 
positions of the clusters with the same redshift distribu- 
tion. Each random catalog has riran points homogeneously 
distributed within the same boundary than the subsam- 
ples. To compute S.cd'r) we have used the estimator 



ecc(r) = 2/ 



7i{r) 



(2) 



where n{r) and nran{f) are the number of cluster-cluster 
and cluster-random pairs separated by a distance r respec- 
tively; / = Nran/ {N — 1) whcrc N and Nran are the total 
number of clusters in the observed sample and random 
catalog respectively. 

We have considered 4 lower limits in X-ray luminos- 
ity which define volume incomplete subsamples of clusters 
(subsamples li to 4i). We have also defined other 4 sub- 
samples by further imposing the restriction of a cut in 
redshift (zcut) in order to build volume complete subsam- 
ples (subsamples Ic to 4c) (see tables 1 and 2). In figure 1 
the horizontal lines define the lower limits of the 4 incom- 
plete subsamples. The regions above the horizontal lines 
and with limiting redshift at the vertical lines define the 4 
complete subsamples. The resulting number N of clusters 
is between 85 and 214 in the incomplete subsamples, and 
between 43 and 72 in the complete subsamples. 



Sample 



N 



Mpc h-^ 



AU 


> 0.02 


248 


-1.92 


21.1 


li 


> 0.27 


214 


-1.89 


22.1 tl-l 


2i 


> 0.54 


168 


-1.80 


23.5 It? 


3i 


> 0.91 


117 


-1.59 


30 1 

•J^--"^ -10.9 


4i 


> 1.38 


85 


-1.75 


^' -16.0 



Table 1. Incomplete Subsamples 

We have fitted the correlation functions obtained with 
power laws of the form £,cc{r) = (r/ro)^. We have es- 
timated the best fitting parameters 7 and tq and their 
associated errors using a a maximum-likelihood estimator 
using a —minimization procedure develo ped by Levem- 



berg and Marquard, (see Press et al. 1987 ). This method 



deals with the errors in each distance bin providing a re- 
liable set of fitting parameters to the correlation function. 
In our calculations we assume Poisson errors c^i n{r) 
in each bin t o estimate the un certainty in the correlation 
function (see Croft et al. 1997 and references therein ). 



Sample 


ix[1044 
h ^erg s 


'] (<) 
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7 


^-o 

Mpc h-^ 


Ic 
2c 
3c 

4c 


>0.27 
>0.54 
>0.91 
>1.38 


0.071 
0.101 
0.130 
0.160 


59 
72 

54 
43 


-1.76 
-1.80 
-1.59 
-1.77 


26.4 +?-^ 

24.6 tti 

30.1 tlo's 

34.7 tin 



Table 2. Complete Subsamples 




Figure 2. Cluster-cluster two-point correlation functions ^cc{r) 
corresponding to the total sample. 
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In figure 2, 3a and 3b arc shown ^cc{i") corresponding 
to the total sample, the incomplete subsamples li-4i, and 
the complete subsamples lc-4c respectively. Error bars in 
^ccir) correspond to Poisson estimates of the uncertainties 
in the number statistics ~ \/n{r). 
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Figure 3a. Cluster-cluster two-point correlation functions 
(,cc{r) corresponding to the different incomplete subsamples an- 
alyzed. 



Figure 4. Ellipses of confidence (1, 2 and 3 a level) corre- 
sponding to the correlation functions shown in figures 2. 




Figure 3b. Cluster-cluster two-point correlation functions 
^cc(r) corresponding to the different complete subsamples an- 
alyzed. 

Estimates of the uncertainties in the power-law best- 
fitting parameters and 7 of the correlation functions 
may be visualized as plots of error contours — Xmz, 
the ro — 7 plane. In figures 4, 5a and 5b we show the 
corresponding error contours of confidence (1, 2 and 3 cr 
level) of the total sample, incomplete subsamples li-4i, and 
complete subsamples lc-4c. 



Figure 5a. Ellipses of confidence (1, 2 and 3 a level) corre- 
sponding to the correlation functions shown in figures 3a. 
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Figure 5b. Ellipses of confidence (f , 2 and 3 a level) corre- 
sponding to the correlation functions shown in figures 3b. 

3. CONCLUSIONS 



a weak increase of the correlation amplitude with in- 
creasing X-ray luminosity which is not statistically sig- 
nificant and suggests a lack of a strong dependence of 
tq on cluster mass. For instance, in our complete sub- 
sample 4i with the highest X-ray luminosity threshold 
Lx > 1.38 10'^'^ h~'^ergs~^ we obtain the highest value 
of correlation length, ro ~ 34.7^jgg Mpc h""'^. Never- 
theless, this value does not differ significantly from tq ~ 
26A~^^'g Mpc h~^ corresponding to subsample Ic with 
> 0.2710-** /i-2er<?s-i. 
There is a well documented evidence for the dependence 
of the correlation length on cluster richness as indicated 
by the relation between tq and the mean inter-cluster sep- 
aration dc in Abell cluster samples. The weak dependence 
of tq with X-ray luminosity threshold as derived from our 
analysis is partial ly related to the broa d relation between 
Lx and richness (Briel & Henry 1994). The relation be- 
tween mass, richness and X-ray luminosity is uncertain 
and is affected by several observational biases and system- 
atics (contamination by projection, departures from hy- 
drostatic equilibrium, etc) as well as astrophysical issues 
(galaxy formation and evolution in clusters, pre-heating of 
the intra-cluster gas, shocks and supernova heating, etc). 
These effects are important for a suitable interpretation 
of the observations given the different mass dependence 
of the cluster correlation length expected in the variety of 
scenarios for structure formation. On the theoretical side 
the situation is also unclear. In hierarchical models of the 
CDM type the dependence of rp o n dc is found eit her very 
weak (Croft fc Efstathiou 1994 ), or moderat e (Bahcal" 



We haw analyrzprl ihp twn-pnint gpatial r-nrrp1atinn fnnr.- 



tion of clusters of galaxies select ed from a, sa.mple of X -ray 



fc Cen 1992), discrepancies that according to Eke et a. 



brightest Abell-type clusters of Ebcling ct al. 1996. For 



the total XBACs sample we find a power-law fit of the 
form ^cc{r) = {r/ro)^ with ro = '2.1.1^2% Mpc h'^ and 
7 = —1.92, values c onsistent with those de rived for R>1 
Abell clusters (see Bahcall & West 1992 and references 
therein). 

In order to provide an insight of the dependence of the 
cluster spatial correlation length ro on mass, we have es- 
timated auto-correlation functions for subsamples of clus- 
ters with different X-ray luminosity thresholds. We find 



1996 may rely on the different cluster identification algo 
rithms. These considerations and the results of our anal- 
ysis suggest that cluster X-ray luminosities may be used 
for a reliable confrontation of models and observations. 
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